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Abstract 

We present a new approach how to calculate the scalar exchange-correlation potentials and 
the vector exchange-correlation potentials from current-carrying ground states of two-dimensional 
quantum dots. From these exchange-correlation potentials we derive exchange-correlation energy 
densities and examine their vorticity (or current) dependence. Compared with parameterizations 
of current-induced effects in literature we find an increased significance of corrections due to para- 
magnetic current densities. 
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I. INTRODUCTION 



Current-spin density functional theory (CSDFT)^ is a powerful tool to calculate ground 
state (GS) properties of interacting electron systems. In contrast to spin-density functional 
theory (SDFT)^ and density-functional theory (DFT)^, CSDFT takes into account the 
coupling of the magnetic field to both, the spin density and the magnetic current density. 
However, CSDFT calculations require the exchange correlation (XC) energy as an input 
which is a functional of the spin densities n a {r) and the paramagnetic current density j p (r). 
In the local vorticity spin-density approximation (LVSDA) the XC functional can be ap- 
proximately calculated using XC energy densities from homogeneous systems which depend 
on the density parameter r s , the polarization £ and the vorticity v (or alternatively the 
local filling factor /). Having investigated the XC energy density for two-dimensional (2D) 
systems with vanishing current density in a previous work^ the focus of this paper is on 
current-induced effects on the XC energy density. The relation^ v = — eBh om /m* allows 
to link the vorticity v = V x to homogeneous systems with a fictitious magnetic field 
Bhom- This property is used by the current- or vorticity- induced corrections of XC energy 
densities in literature. They are parameterized as function of Bh om or the local filling factor 
/ = 27rn?B- The local filling factor defines the number of occupied Landau levels of a 2D 
system in a magnetic field B. l B = y / l/(e|B hom |) is the magnetic length corresponding to 
Bhom (h = 1 throughout this paper). 

CSDFT calculations for 2D systems (for quantum dots (QDs) see e.g. Refs. PJflflcl]) 
have made use of current-density dependent XC energy densities based on interpolations 
provided by Levesque, Weis, and MacDonald (LWM)^ or Fano and Ortolani (FO)^. Both 
parameterizations are based on the assumption that the coupling of the magnetic field to the 
currents is only relevant for high B. For these strong magnetic fields the considered system 
is in the lowest Landau level as the inter-Landau level spacings increases with growing B 
and the Zeeman energy and exchange effects favor spin polarization. 

For this regime Laughlin^ postulated the relation 

e xc (f = m- 1 ) = --= dx[g f (xl B )-l] (1) 
Ibv 2m Jo 

for filling factors / = m _1 < 1 (m is an odd number). For discrete filling factors f — \i \ 
Levesque, Weis, and MacDonald (LWM)^ calculate the XC energy density by using the pair 
correlation function g/(r) of a classical plasma with a density n = //(27r/ B ) and statistical 
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methods.— Interpolation delivers a continuous expression for the dependence of the XC 
energy density on the filling factor / which is valid for / < 1 

^xc M (^ /)/Ry = -0-782 133 — (l - 0.211 / a74 + 0.012 / L7 ) . (2) 

Fano and Ortolani (FO) 1 ^ pursue an alternative concept to parameterize current correc- 
tions for exc'- They follow the particle- hole symmetry for electrons in the lowest Landau 
level at large magnetic fields and establish an interpolation which exactly reproduces the 
well-known analytical result for filling factor / = 1 (see Ref. Q|) 
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(3) 



In spite of concerns that the exact XC energy density might be irregular and not 
diflerentiable,— the parameterizations (J2J) or © are used in a number of paper a 8 i 9 i 10 i 17 as- 
suming that they reproduce an important part of current-induced XC effects. However, 
the LWM and the FO parameterization are only correct in the high density limit, i.e. for 
r s -> O.ii 

To obtain XC energy densities Exc{ r s,£,, /) these current-dependent corrections are com- 
bined with XC energy densities £xc^ LWM ^ F ° ( r s, from systems without current- density by 
an interpolation established by Rasolt and Perroij 1 ^ 

LWM/FO, n fp TC/AMBG/ISI, n 

exc(r s , Z, f) = ^ ^di±L^ . (4) 

For large filling factors / (i.e. for small magnetic fields Bh om ) the Pade approximation yields 
the well-known result of local spin density approximation (LSDA), £xc / ' AMGB//ISI ( r s; 0- For 
this quantity the interpolation by Tanatar and Ceperley (TC)^ is the most frequently used 
approach. However, recently new parameterizations were proposed by Attaccalite, Moroni, 
Gori-Giorgi, and Bachelet (AMGB)^ and Seidl (ISI=interaction strength interpolation)^. 
In the limit of high effective magnetic fields or / < 1 Eq. (plj takes into account corrections 



due to paramagnetic currents. Usually the exponent is chosen to be p = 4 (see Ref. lsl]). 
However, results seem to be almost insensitive with respect a variation of p . 17 ' 18 The usage 
of the interpolation (j3J) is justified by an improvement of CSDFT results compared with 
benchmark calculations from ED and Quantum-Monte-Carlo.^ 7 - 



It is also reported that results do not significantly depend on whether the LWM or the 
FO parameterization is applied.— Even working with corrections which take into account 
discontinuities in the derivative of the XC energy density with respect to / in the fractional 
quantum Hall regime (see Ref. (2^]) or an additional dependence on the polarization (see 
Refs. QQ) 

do not further improve the CSDFT results.— 
The present paper is a natural extension of our previous publication^ where we presented 
2D XC energy densities for systems without paramagnetic current density which were ex- 
tracted from exact GS densities and energies provided by exact diagonalization (ED) 2 4. Here 
we want to focus on the current-induced effects on XC energy densities. Analogous to Ref. 
[3] the extraction of XC energy densities is performed in two steps: After briefly introducing 
CSDFT in Sec. II and the properties of the QD Hamiltonian in Sec. Ill we present the 
first step of the method, i.e. the calculation of scalar and vector XC potentials from exact 
QD GSs together with results (Sec. IV). In the following Section we introduce the second 
step which comprises the concept how to obtain XC energy densities from scalar and vector 
XC potentials. The results are compared with the current-dependent XC energy densities 
discussed above. Sec. VI gives a summary of the paper. 



II. CURRENT-SPIN-DENSITY FUNCTIONAL THEORY AND LOCAL VORTIC- 
ITY (SPIN-)DENSITY APPROXIMATION 

In this Section we sketch the basics of CSDFT and LVSDA. The DFT formalism was 
originally established by Hohenberg, Kohn, and Sham^ and generalized to spin-polarized 
systems^ and current-carrying systems^ by including the coupling of the polarization to an 
applied magnetic field and to the paramagnetic current density. Accordingly, the Hohenberg- 
Kohn (HK) theorem has to be modified with respect to the spin degrees of freedom^ and 
to the paramagnetic current density^. For the most general case, it states that two dif- 
ferent non-degenerate GS wavefunctions and \^') always yield different combinations 
(n a (r), j p (r)) 7^ (n^(r), j' (r)) of spin densities and paramagnetic current densities. This is 
sufficient to establish a functional of the total energy with the usual functional properties 

Ev 0> An a ,j p ] = FhkKJ p ] + E / drV;(rK(r)+e J dr j p (r)- A(r) + ^ ]T |dr n CT (r)A 2 (r) 
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and the universal HK functional 

^hkKJp] = (^K,j p ]|T + iy|^K,j p ]). (6) 

However, in contrast to the original HK theorem^ the uniqueness of the potentials is 
lost in CSDFT, because it is possible that two different combinations of spin-dependent 
potentials and vector potential yield the same GS wavefunction. Thus, the one-to-one map 
between GS densities and XC potentials is destroyed. The consequence is that there might be 
discontinuities in the functional derivatives of the HK functional or the XC energy functional 
with respect to the densities. 26 Capelle and Vignale 2 ^ conclude that this problem also affects 
the extraction of exact XC potentials from given GS densities. However, in the case of axially 
symmetric, parabolic quantum dots it is possible to restore the uniqueness as we will see in 
Section H3 

For practical purposes the variational scheme has to be mapped to the Kohn-Sham (KS) 
system, i.e. an effective single-particle system with the same GS densities as the interacting 
system. The spin-degree of freedom is considered in the KS equations 3 by assuming the 
total spin S z in z-direction to be a good quantum number 

-L {Z + e [A (r) + A xc (r)]} 2 + ^ {A 2 (r) - [A(r) + A xc (r)] 2 } 



+ V a (r) + V H (r) + Vi 



xc,, 



ip j; a{r) = £ j!(T ip j;a {r) (7) 



with the spin a = ± in z-direction and the KS energies Ej )(T in increasing order. For a system 
containing N particles the spin densities and the paramagnetic current density are given by 

M r ) = Et^I^MI 2 ( 8 ) 

3 

JpOO = hVWV^r) - (v^* CT (r)) ^(r)] (9) 

with 7 Ji(T being occupation numbers of the KS levels in the GS (J2j lj,a — N a and A^ + + = 
N). Then the GS density and polarization are 

n(r) = n+(r) + n_(r) (10) 

fw = " +(r) 7"- (r) (ID 

n(r) 



The XC potentials 



VxcAKM = (12) 



eAxc, CT (K,j p ],r) = — — — (13) 

°Jp[ r ) 

are defined as functional derivatives of the XC energy functional 

#xcK, jp] = FhkK, j p ] - y dr y dr'^-^i - T s [n a , jj. (14) 

(Ts[n CT , j p ] denotes the kinetic energy functional of the KS system.) The total GS energy S 
of the interacting system can be calculated from 

= S^>^>~54^/ dr / dr '"yZ~^~5I y"drV r xc,a(K,jp],r)n CT (r) 
-e / dr jp(r) • A xc (r) + E xc [ "(j) JpJ • 

(15) 

Concerning the XC potentials we apply the LVSDA 

^xcK,v]~ J drn(r)e X c(n a (r), v(r)). (16) 

Using the vorticity v(r) = V x instead of the paramagnetic current density is 
a consequence of the gauge invariance of CSDFT and appropriate for local density 
approximations. 1,2 The most important parameterizations of the XC energy density used 
in 2D calculations were introduced in Sec. I. 



III. QUANTUM DOT HAMILTON! AN AND GROUND STATE DENSITIES 

We consider a two-dimensional QD with an axially symmetric parabolic confinement 
potential of strength u in a magnetic field B = (0, 0, B). The Hamiltonian for N particles 
in real-space representation (with r = (x,y),p = {p x ,Py), angular momentum in ^-direction 
l z = xp y — yp x = and vector potential A(r) = y, x)) reads 

* = % to + + + i|; do 

or 

# = V f-^- + -m*u 2 h r] + ^0 + -f f- -. (18) 

3 eA2m* 2 h > 2*) 2^ 47ree |r i -r fc | 1 ; 

Here e is the electron charge, m* = 0.067 m e and e = 12.4 are the effective mass and the 
screening constant of the host semiconductor (assumed to be GaAs), u c = ^4 is the cyclotron 



frequency, and uih = J + cu 2 /4 is the hybrid frequency. 



The reference densities and energies for the investigation of XC energy densities are 
calculated by ED techniques which provide results of high accuracy— Note that the (spin- 
density and the current density of all eigenstates of the angular momentum operator are 
functions of radius r but not of the angle (p. Moreover, the current density is always of 
the form j p (r) = j p ^(r) e 9 whereas the radial component vanishes. Therefore, the relevant 
quantities provided by ED are the spin-densities n a {r) and the azimuthal component j P)l p(r) 
of the paramagnetic current density (or the ^-component v of the vorticity) of the GS and 
its energy E . 

The many-particle Hamiltonian also shows an interesting scaling property^ which plays 
an important role in the context of the HK theorem. It is revealed if we rewrite (I18|) as a 
function of ujq, uj c , and iV and separate the angular momentum contribution. This concept 
is similar to the modification of the Hamiltonian of a (natural) atom in Ref . 2(| : 

H(lo 0i uj c ,N) = V ( —p 2 ; + -m* ( ul + 4] ri\ + - E' f- r + E ~k 



UJr 



= H(^u 2 +u*/i,0,N) + ^L. (19) 
We identify the first part of (pm)) as a Hamiltonian of a QD in zero magnetic field with the 



confinement frequency y Uq + The sum over all single-particle angular momenta in 

the last term yields the total angular momentum. As a result we can map the spectra and 
wave functions of QDs in the magnetic field to QDs in zero magnetic field. 
The KS Hamiltonian of axially symmetric 2D QDs reads 



2m* 



I f r d_\ I 2 



r dr V dr I r 2 dip 2 



+ \m^y + V H {r) + V X cAK, v], r) (20) 



Obviously, the KS Hamiltonian shows the same scaling properties as the original Hamiltonian 



H s>a (u; , lu c , N) = H s ^u 2 + u; 2 /4, 0, N) + y L. (21) 

Therefore, the XC effects do not depend on the strength of the external magnetic field. We 
can restrict our studies to systems without magnetic field. However, without magnetic field 
we can apply the original HK theorem to the Hamiltonian (fTTj) what guarantees the unique- 
ness of the wavefunction and of the external potential. For the KS system the uniqueness 
of the (XC) potentials in quantum dots is explicitly investigated in the following Section. 



7 



IV. SCALAR AND VECTOR XC POTENTIALS 



A review of the literature on the numerical inversion of KS equations, i.e. the calculation 
of XC potentials from GS densities can be found in Ref. []]]. However, in this paper we 
will focus on an unpolarized and on a fully polarized system for which we can analytically 
calculate the scalar XC potential and the XC vector potential. Before we present our results 
for Vxc and Axe we will discuss the problem of uniqueness of XC potentials in CSDFT in 
the next paragraph. 



A. Uniqueness of XC (vector) potentials 

Because of Eqs. (fTTJj) and (|2*TJl we can assume the single-particle states \(pj >a ) (j is a com- 
bination of radial and angular momentum quantum number) of the KS Slater determinant 

1^) = l^u>-->^a>l^i,T>---I^Jv T >t) ( 22 ) 

to be eigenf unctions of the KS Hamiltonian ()20|) at B = 



Hks = 



2m* 



I f r _d_\ I 2 



r dr \ dr I r 2 dip 2 



1 * 2 2 , d e A *cA r ) 

+ -m uj r + 



2 ' dip im* r 



+Vu(r)+J2y^Ca(r)\cT)(a\. (23) 

a 

with the potentials (Vxca, Akc i¥ >) 

H KS \ipj,o) = £j>|^>)- (24) 

Due to the axial symmetry of the system the wavefunctions \<Pj >a ) are a l so eigenfunctions of 
the angular momentum operator l\(pj >a ) = lj,a\<Pj,<r) ■ 

Now let us assume that there is a KS Hamiltonian if^s with the potentials 
(Fxcd ^xc,</>) 7^ (Vkco- + c a, Axc,tp) so that the KS single-particle wavefunctions are also 
eigenfunctions of 

#ksI^> = 4,>^>}- ( 25 ) 

Subtracting Eq. (f2"2j) from Eq. (f2~5|) delivers the constraint 

AVxcv(r) + k^ ^ /^ - 4. + ^ = ( 26 ) 



S 



for the potential differences AVxc,*(r) = V XCa {r) - Vxc,*(r) and AA X c, v ( r ) = A 'xc, v ( r ) ~ 
Axc,ip( r )- Eq. (|2l)Jl has to be fulfilled for all wavef unctions \<fj,a)- 

If at least two states with different angular momenta lj >a ^ ly^ are occupied we can 
conclude from 



that the scalar potentials agree up to a constant AVxc,o-( r ) = const. Thus, AAxc )ip {r) = 
const' • r immediately follows. As a consequence of finite XC vector potentials at B = we 
derive AAxc,ip{f) = 0. If there are also occupied states in the other spin direction we obtain 
AVxc-a{r) = const. This result is in contradiction to (V XC ^,A' XC ^) ^ (Vxca + c a ,A xc , v ) 
and the uniqueness of XC potentials is proven. 

The assumption of at least two occupied states with different angular momenta lj >a 7^ lj^ a 
in one spin direction is correct for the configurations considered in this paper. All other cases 
are more or less pathological: 

• For iVf = iVj_ = 1 uniqueness cannot be proven except for the the case = l±i = 
(i.e. conventional DFT without paramagnetic current densities and vanishing Axe)- 

• In the case of iV > 2, /ij = ... = ijv.j, and Zij = ... = In^, uniqueness cannot be 
guaranteed either. However, such a configuration should not be a GS in particular for 
larger particle numbers. 

In general the XC potentials are overdetermined by ()2fij) . Therefore, the crucial aspect is 
the existence of appropriate XC potentials but not their uniqueness. 

B. Results for XC (vector) potentials 

From the discussion of the uniqueness of XC potentials presented above immediately 
follows that we can analytically calculate the scalar and the vector potentials if we are able 
to derive the exact KS wavefunctions from (n CT (r), j p (r)). This can be accomplished for 
special systems. 

The simplest conceivable configuration is the fully polarized two-electron v = 1-quantum- 
Hall-droplet (QHD) with quantum numbers L = — 1, S — 1, and S z = — 1. In particular the 




(27) 
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confinement energy uj^ = 2/3Ry allows for a completely analytical solution.— (For the rest 
of this subsection we denote energies in Ry and lengths in effective Bohr radii ao.) According 
to Taut— the separation of center-of-mass and relative motion delivers the GS wavefunction 

#(ri,ri) = v9 C M(|ri + r 2 |/2)v9 rc i(|r 1 - r 2 |) 

= C CM e- R2/3 C rcl r (1 + r/3) e"^ 12 (28) 

(R = (ri + r 2 )/2) and r = ri — r 2 ) and GS energy 8/3. The constants are Cqm = and 



C rel = 1/^2tt [9V6V + 42J. 

Using these results we can analytically calculate the GS density n(r), the azimuthal com- 
ponent of the paramagnetic current density j PtV ,(r), and the (^-component of the) vorticity 

v{r) = -— r — — — . (29) 
r or \ n(r) J 

Now we are interested in the effective potentials which exactly reproduce the GS density 
combination (n(r) , j p ^(r)) . In this context we make use of the scaling relations (jl9)l and 
(l2*Tj) and restrict ourselves to the case B = 0. 

Labeling states and occupation numbers by three indices which characterize the radial 
quantum number, the angular momentum and the spin state, the GS of the (spin-down) 
v = 1-QHD in the KS system is given by 70,0,-1 = 7o,-i,-i = 1- Then the relations between 
KS wavefunctions (fijj^(r) = -i=e im¥ 'i?j ) i(r) and density (Fig. da)) 



"W^^^Wfe^W (30) 

and paramagnetic current density (Fig. [TJb)) 

j P A r ) = J2 ^jM^hA^hiA*) ( 31 ) 

allow for the calculation of the vorticity (Fig. ^c)) and of the exact KS orbitals 



iV(r) = y/^/Ry)rj PtV (r) (32) 

and 



Ro, (r) = y/2Tm(r) - ti% tl (r). (33) 
The radial KS wavefunctions i? ,o( r ) un d -Ro,i( r ) are depicted in Fig. Hfd). 
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Consequently, we can calculate the effective scalar single-particle potential up to a con- 
stant £0,0,-1 (which is a KS energy) by inverting the KS Schrodinger equation for 0o,o,-i(r) 
(see Fig. [He)) 

)& (r&)] i?o,o(r) 



K ff (r) = r 2 /9 + V u (r) + V xc (r) - £o ,o,-i 



^o,o(r) 



(34) 



In particular the KS equation for states with I = only involves the scalar potential, but not 
the XC vector potential. As for given density n(r) the Hartree-potential Va(r) (Fig. ^f)) is 
also known we can derive the exact scalar XC potential (Fig. HJg)) 



V xc {r) =£ ,o,-i -r 2 /9- V H (r) + 



Ro,o(r) 



(35) 



up to the gauge constant £0,0,-1 • 

Finally we can find an approach to the XC vector potential (Fig. ^h)) by exploiting the 
KS equation for o _ 1 _ 1 (r). After solving the KS equation for Axc,<p{i") an d substituting the 
scalar XC potential we obtain 



1 r 
Ry2 I r 2 



Ro,i(r) 



Ro,o{r) 



+ £0,0,-1 — ^0,-1,-1 



(36) 

The singularity in the term 1/r 2 is eliminated by the second one as both terms together 
represent the kinetic energy density of the state 0o,-i,-i( r ) which is finite at r = 0. 

Now let us focus on the asymptotic of all relevant quantities. In the limit of large radii 
the asymptotic expansions of the densities and the vorticity are 



limn(r)^|c7§ M ^ el r 4 e- 2 /3 
unj^r) ^ 3(1/Ry) 



lim v(r) — ► Ry (-24/r 4 + 54/r 
Thus, the results for the radial wavefunctions i?o.o( r ) an d -Ro,i( r ) are 



lim # ,o(r) -> 2V27cC C uC vcl re~ r2 / 6 

r—>oo 



lim R i(r) 



27rCcMCrel 2 -r 2 /6 

= r e ' . 

\/3 



The effective potential is parabolic in the limit of large r 



lim VLffO 



1 2 4 3 3 

-r 1 1 . 

9 3 2r 4r 2 



(37) 

(38) 
(39) 

(40) 
(41) 

(42) 
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With the Hartree potential converging to zero 

lim Vkfr) -> - (43) 

we assume that the XC potential vanishes in the limit r — > oo. 28 This is achieved by an 
gauge constant £0,0,-1 = 4/3 

limV xc (r)^-|- + A. (44) 

The asymptotic of the XC vector potential gives a constraint for £0,-1,-1 = 2 if we do not 
want Axctpir) to diverge in the limit of large radii. (This is a reasonable requirement as we 
work at B = 0. A contribution Axc,ip( r ) — %;T would correspond to an external magnetic 
field B'= (0,0,5').) 

lim A xc » -> J- (-- - pi . (45) 

Thus, the XC potentials {Vxc{ r ), ^4xc,^( r )) are uniquely determined. 

Besides this fully analytic approach it is also possible to perform the calculation of XC 
potentials starting with GS densities from ED. Using numerical data for the densities allows 
us to study the dependence of the XC vector potential on the strength of the external 
confinement potential. In Fig. |21 we show the scalar XC potentials (up to a constant) and 
the XC vector potentials for ujq in the range between 3 meV and 100 meV. In the context of 
the XC vector potentials we find a pronounced sensitivity to the precision of the densities 
at the edge of the dot. Because the Axc#>( r )) is calculated as a difference of two (large) 
effective potentials, numerical errors are growing at the edge of the dot although 16 Landau 
levels were used in ED. However, the results from the bulk part of the QD are reliable and 
show a continuous dependence on the strength of u . 

The second system for which we study scalar and vector XC potentials is the unpolarized 
four-electron system [y = 2-QHD) with S = 0, L = —2. In this case all GS densities are 
provided by ED using 335259 Slater determinants (see Fig.|3]for GS (paramagnetic current) 
densities at ujq = 3meV). Concerning the mathematical structure of the KS equations it is 
similar to the v = 1 two-electron droplet: the v = 2 four-electron droplet can be regarded as 
a combination of two v = 1 two-electron droplets with different spin polarizations. Therefore 
the inversion of KS equations is completely analogous to Eqs. (|3U |) -(|3l) ]) if we use half of the 
GS density and paramagnetic density of the four-electron v = 2-QHD. 
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Fig. |3] shows the results for scalar XC potentials (up to a constant) and vector XC poten- 
tials as a function of the strength of the confinement potential. Concerning the precision of 
the vector potentials we can repeat the same arguments as for the two-electron system, i.e. 
only the values from the bulk part of the QD are reliable. Similar to the v = 1-droplet we 
find again a continuous dependence on u . Another aspect is the size of the amplitude of the 
XC vector potential which is even larger than for the two electron system (see Fig- EJ) - This 
represents a strong hint that current-induced corrections are for the v = 2-QHD at least as 
important as for the v = 1-QHD. This aspect is not considered in calculations relying on an 
interpolation of the type of Eq. (J3J) which neglects the influence of higher Landau levels. 

V. VORTICITY-DEPENDENT XC ENERGY DENSITIES 

In this Section we will deal with the vorticity dependence of XC energy densities extracted 
from the XC scalar and vector potentials above. The two parts comprise the methodology 
and the results for the four-electron v — 2 system. 

A. Calculation of XC energy densities 

The steps for extracting XC energy densities from XC potentials are analogous to Ref. 
0] but generalized with respect to the vorticity dependence. Eqs. (JT2J) and (fE^) provide the 
general relation between XC functional and XC potentials. To construct a relation between 
XC potentials and XC energy densities we apply the LVSDA (fTBj) and obtain 




(46) 



VxcA r ) = £ xc(r) + n(r) 



dexc 



(r) - eA X cA r ) 



Jp,<p\ r 



n(r) 



(47) 



We solve for ^gf (r) and §^(r) 




(49) 



(48) 



considering the (physical) boundary condition 



r— >oo 



lim n(r 



de X c 
dv 



(r) = 0. 



(50) 
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The results are substituted in the derivative of the XC energy density as a function of the 
radius 



d£xc(V) _ y> 9e X c dn a de X c dv 
dr Q dn a dr dv dr 

After some transformations we arrive at a linear differential equation (DEQ) 



dr 



+ 



dr 



(rjexcvr) = I[r) 



(51) 



(52) 



with the inhomogeneity 

1 



I r) 



^ , E (W) + *) d ~^(r) + eA xcA r) ^ dh " n (r) 



1 dv , . 
n(r) dr 



dr"n(r")eA xc Jr") 



(53) 



which contains all the information about the dependence on the spin densities and the 
vorticity. In Eq. (JHHj) we take into account that the gauge constants c CT of the scalar potentials 
are not known from the previous step. They will be calculated later. The solution of 
the homogeneous part of DEQ (J52J) is given by ^xc 1 ! 7 ") = ®/ n { r ), a special solution can 
be calculated using the ansatz £xc°( r ) = /3{r)/n(r). The function (3{r) follows from an 
elementary DEQ 

(54) 



M( r ) = n (r)I(r), 
dr 



whose solution is 



Thus the general solution of DEQ (J52*j) is 



^xc(r)=4°c»+4 P c C W 



/3(r) = / dr'n{r')I(r') -/3(0). 



(55) 



a -(3(0) 
n(r) ' n(r) 



+ 



dr' n(r')I(r'). 



(56) 



After substituting the inhomogeneity the result reads 



a -/3(0)- n a (0)J2aCa . M r ) 

£xcl r J — ~77\ r 2^ c <y 



n(r) 



n(r) ' ^)/ dr '{? Fx ^ (r,) ^ (r,) 



9 log 7T (9f 

+eAxc, v (r , )J P ^(r , )^^(r , ) + ^(r ) 



dr // n(r // )eA x c, ¥ >(0 
With the modulus of the XC energy being finite in LVSDA 



(57) 



Exc[n a ,v]\ « y dr n(r)£ XC (r) = / dr n(r)|e X c(^)| < oo 



(58) 
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we choose a — (3(0) — n a (0) J2a c a — thus avoiding any divergent contributions. Conse- 
quently the analytical solution of (J52j) satisfying the physical boundary conditions is 



£xc(r) 



„ n(r) n 



l V) /d^EWO^M +eA^(r>)j P , 9 (r') 



, s d\ogn, 

r 



9r' 



&r"n(r")eA^(r") 



(59) 



The last step is the calculation of the gauge constants c a of the scalar potentials. One 
condition which has not been used yet is the agreement of the DFT GS energy and the 
exact GS energy. If we eliminate Vxc,o-( r ) i n Eq. (JT5j) and write it in a modified form 

E GZ = Et,>/ dr^V(r) (-L ( p + e A(r)) 2 + K(r)) ^(r)+± / drV H (r)n(r)+E xc [n a ,v] 

(60) 

we can calculate the (exact!) XC energy Exc{ n a, v ] f° r GS densities: the KS wavefunctions 
are known from the selfconsistent solution of the KS equations. Thus the expectation values 
/ dr(p* a (r) (p + eA(r)) 2 + V a (rfj </^>( r ) of the non-interacting system can be calcu- 
lated. The Coulomb energy and the GS-energy are uniquely determined by the interacting 
system. On the other hand the XC energy in LVSDA is given by 

E X c[n a ,y] ~ J drn(r)e X c(r) 

= 'EcN.+ Jdr jdr4^ V xcAr')^(r') + eA xc , v (r') x 

a r=0 n ^ °" 



. , ,,d\ogn A 



dv 



dr"n(r")eAxcA r ") 



(61) 



dr' dr' 

what makes J2 a c aN a accessible. In the case of unpolarized systems c := c-f = cj_, and 
for full polarization the constant c CT of the unoccupied spin direction is irrelevant (as the 
corresponding N a = 0). Consequently, the results for the XC energy density are unique. 
For partially polarized systems uniqueness of the results can be achieved by additionally 
demanding asymptotic agreement of Vxc,o-( r ) for finite systems^ 

lim VxcM = lim V X c,i(r). (62) 

r—*oo " >* v 

As a result of the inversion of the LVSDA formalism we obtain the XC energy density 
as a function of the radius. After eliminating the radius using the spin densities n a (r) and 
vorticity v(r) or alternatively r s (r), £(r), and v(r) we arrive at the representation of the XC 
energy density as a function of (r s , £, v). 
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B. Results for vorticity-dependent XC energy densities 

In this Section we summarize the numerical results for the extracted XC energy densities. 
Because of the high contribution of self-interaction in the two-electron system only the 
four-electron v = 2-QHD can be used for the derivation of vorticity-dependent XC energy 
densities. After solving the DEQ (|52|) for different confinement energies ujq we get a set 
of 6xc f° r different r s , v, and £ = (no polarization). The results for different vorticities 
are shown in Fig. El as a function of r s together with the pure TCUS parameterization, TC 
combined with LW M'W 9 , and TC combined with F Q 12 ' 18 ' 19 . Whereas the LWM and 
FO parameterizations do not deviate very much from the pure TC data which does not 
consider any current- induced effects, the extracted XC energy densities exhibit a pronounced 
dependence on the vorticity. In particular, the curve of extracted exc starts for small r s 
below the reference curves, crosses it and finally lies above it. Moreover, the crossing wanders 
to larger r s if we decrease the vorticity (see Fig. Efa)-(f)). 

Although the interpretation of these large deviations from the standard parameterizations 
is difficult as the data set is very small our results give rise to the question if the vorticity 
dependence of XC energy densities can be treated in the frame of LVSDA. There might 
be also non-local effects from paramagnetic currents which cannot be considered in this 
approach. However, in any case it is necessary to reexamine the vorticity-dependence of exc 
with focus on Landau level mixing and arbitrary densities. 

VI. SUMMARY 

In this paper we presented a new approach to study scalar XC potentials and XC vector 
potentials of current- carrying ground states in axially symmetric, parabolic quantum dots. 
Starting from exact ground states densities and paramagnetic current densities we analyt- 
ically calculated the XC potentials of the two-electron v = 1-droplet and the four-electron 
v = 2-droplet and studied them as a function of the strength of the confinement potential. 
In both systems the XC vector potentials are of the same order of magnitude which shows 
the importance of this current-induced effect even for higher filling factors. This result is 
not taken into account in the interpolation introduced by Rasolt and Perrot^. 

The second step involved the extraction of vorticity-dependent XC energy densities. 
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The effect of the vorticity is far more pronounced than predicted by Rasolt and Perrotf^, 
Levesque, Weis, and MacDonalcUi and Fano and Ortolani 12 . A possible reason for this be- 
havior is the assumption that the influence of current densities is only relevant for strong 
magnetic fields when the local filling factor / is less than one. This is the restriction which 
enters the derivations of the parameterizations in literature. However, our results point 
out that a new approach to current- or vorticity-dependent XC energy densities is required 
which also comprise the effects of Landau level mixing. 
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FIG. 1: Analytical results for the two-electron v = 1-QHD with confinement energy 2/3 Ry. 
Figs, (a)-(h) depict the GS density, the paramagnetic GS density, the vorticity, the radial KS 
wavefunctions, the effective scalar potential, the Hartree potential, the scalar XC potential, and 
the XC vector potential. 

FIG. 2: Scalar XC potentials (a) and (c) and XC vector potentials (b) and (d) of the two-electron 
v = 1-QHD as a function of the strength of the confinement potential. In (a) and (b) the results 
for luq = 3 meV - 10 meV (stepsize 1 meV) are shown, in (c) and (d) results for ujq = 10 meV - 
100 meV (stepsize lOmeV). The arrows indicate curves corresponding to increasing confinement 
potentials. 



FIG. 3: GS density (a) and GS paramagnetic current density (b) of a four-electron v = 2-QHD 
with confinement energy 3meV. 
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FIG. 4: Scalar XC potentials (a) and (c) and XC vector potentials (b) and (d) of the four-electron 
v = 2-QHD as a function of the strength of the confinement potential. In (a) and (b) the results 
for ujq = 3 meV - 10 meV (stepsize 1 meV) are shown, in (c) and (d) results for uoq = 10 meV - 
100 meV (stepsize lOmeV). The arrows indicate curves corresponding to increasing confinement 
potentials. 

FIG. 5: Vorticity-dependent XC energy densities extracted from the four-electron v = 2-QHD. 
(a)-(f) show the XC energy densities as a function of r s for different vorticities. 
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Fig. 2, Wensauer, PRB 
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